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We study the excited states of the Q-balls by performing stationary perturbation on the spherical 
Q-balls. We find the exact solution of the stationary perturbation of the global Q-ball with thin wall 
approximation. For local Q-balls we solve the equations of motion for the perturbative part approx- 
O |' imately by using expansion about the coupling constant. Furthermore we comment on the magnetic 

field generated by the excited Q-balls during the phase transition precipitated by solitosyhthesis and 
give an implication into cosmology. 

> : 

■ I. INTRODUCTION 

The existence of the coherent magnetic field over various astrophysical scales has been established up to the present 
0. To explain the observed magnetic field for galaxies, clusters and inter-clusters, we need a generation mechanism 
of the seed of primordial magnetic fields || . This is because magnetohydrodynamics shows that the magnetic fields 
cannot exist unless there is non-zero field strength as a initial condition. 

It is known that if the seed of magnetic fields with 10~ 19 Gauss exists over the present comoving scale of the proto 
galaxy, ~ lOOkpc, these fields might be amplified to the observed galactic magnetic fields, 10~ 6 Gauss, by the dynamo 
mechanism || [| 

Recently, several generation mechanism of the magnetic field in the electroweak phase transition have been proposed. 
One of them is based on the thermal fluctuations || . However, this mechanism cannot provide enough coherent scale 
of the magnetic field, and the net magnetic field will be too small to explain the present observation. Another 
mechanism is founded on the bubble nucleation Q. This mechanism is able to supply the large coherent length of 
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the magnetic field and it can provide the enough field strength for the onset condition of the dynamo mechanism. In 
this case, however, one needs to assume that the phase transition is strongly first order one and completed by the 
expansion of the critical bubbles. Although several attempt to clarify the detail of the phase transition have been 
performed, the order of the phase transition has not been clear at this stage m. 

In this paper, we considered the magnetogenesis due to excited states of Q-balls during the phase transition, 
suggested by one of the present authors ||. Q-ball is non-topological soliton solution of the complex scalar field 
arising in the theory with an unbroken continuous symmetry j9l JlOj . The generation of Q-ball jnj and the possibility 
of the phase transition precipitated by solitosynthesis Jl2| have been actively investigated. These studies are 
based on the ground state of Q-ball. If we consider the excited states of Q-ball, however, it is conceivable that Q-ball 
has angular momentum. Then, it could also have the magnetic moment if it is the gauged Q-ball [fl4j| . 

We re-analyze the excited state of the Q-ball by the stationary perturbation analysis, not by the surface wave 
investigated by Coleman ||. We should note that the stationary perturbation is not general. However, it is still worth 
to investigate them because one may expect that the thermal distribution is decided by stationary states. We found 
that the £ = 1 mode which does not exist as in the case of surface wave. The result might have some important effects 
on cosmology, since the contribution of the I — 1 mode is larger than the model t — 2 in general. Furthermore, taking 
account of the above fact and assuming copious production of gauged Q-balls, we estimate the magnetic field of the 
excited Q-balls. 

The rest of the present paper is organized as follows. In Sec. II, and III, we analyze the stationary perturbation on 
the global Q-balls and local Q-balls of the ground states and estimate the conserved quantities such as the angular 
momentum and the magnetic field etc. In the subsections of the above two sections, we review the feature of the 
Q-balls in ground state, and we also give a sketch of the phase transition due to gauged Q-balls. Finally, we give an 
implication into cosmology in Sec. IV. 



II. GLOBAL Q-BALLS 



A. Review of Q-Balls in the Ground State 

In this subsection we review briefly global Q-balls in the ground state. For simplicity, we consider only the complex 
scalar field with a global U(l) symmetry. According to the Coleman's study, the Q-ball has stationary and spherical 
configuration, 

h = Mr)e iut . (2.1) 

For the existence of the Q-ball solution, we consider here the potential which has two minima, that is to say, true 
vacuum and false vacuum. In the thin wall approximation, the spatial configuration of the scalar field, y>o{r), is 
characterized by the radius of Q-ball, R, the wall width, <5, and the field values at the center of configuration, er+. 
Here we set the field value at the minima: 0, cr+, and the top of the potential barrier: er_. Then the total energy can 
be written by 



2 J (p 2 d 3 x 2 



where U± := U(a±) and Qq is the total conserved charge: Qq = J d 3 xipQ. 

The most favourable configuration is determined by the variational principle to minimize the total energy. From 
dE /dS\s=s, — 0, the wall width has the expression 

y/2 JTJZ 



Thus, the total energy becomes 

£0* = E Q (S*) = -^-JrV + 4\/27r<7 + ^/TjZR 2 + ^-U+R 3 . (2.4) 
If the potential has two degenerate minima, U+ = U(0) = 0, the total energy becomes 
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From dE /dR\n = ji_ M — 0, one obtains the radius at which the energy takes minimum value, 



V64\/27r 



D* 



2/5 



Then, the total energy becomes 



E - E *{R ) - — g J 



,/ ^3/5^1/10- 
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For simplicity, we write i?S as R hereafter. 



(2.5) 



(2.6) 



(2.7) 



B. Stationary Perturbation on Q-Balls 

Let us consider the stationary perturbation on the Q-balls. The perturbation was analyzed in Ref. [jl0| from the 
viewpoint of the stability. In this paper, we give explicit solutions of the stationary perturbation. Although we 
will not give any implication to cosmology in this calculation on the global Q-balls, the analysis give a educational 
examination for the extension into the local Q-ball in which we can give an implication into cosmology. 

The Lagrangian of the complex scalar fields is 



L 



d A x 



\\<i>\~lm 2 -uM) 



where t^(|0|) is the potential. The energy momentum tensor becomes 



T^v = -V M fV,, 



U{\d>\) 



We add the stationary perturbation on the Q-ball in the ground state as follows, 

= 00 + 0! = (y> ( r ) + V4^^ m (r)r £m (6»,(^))e^ t , 

where I > 1 and I > \m\. Then the Lagrangian becomes 
L = L Q + Li 



= I d 3 x 



\»W - i(V^o) 2 - Ufa) 



r / d 3 x 



V dr 



- V e (r)(ip er , 



where 



V t {r) = U"(<po) -u 2 + , ■ 
One can see easily from the eq. (2.11) that the perturbation ipi m satisfies 

d 2 ip£ rn 2 dtp i 



dr 2 r dr 



+ Vi(r)tf lm = 0. 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



The above equation is similar to the Schrodinger equation with the potential Vi(r) in the ground state. In the thin 
wall approximation, one can set 

U"(Mr)) = ^(r -R) + li 2 0(R - r), (2.14) 
where /Zq = U"(0) and — U"(cr + ), and then the eq. (2.13) can be solved exactly. For r < R, the solution is 



3 



<Ptm{r) = vf(r) = A' e j e (x) 

\r dry \ r / 

where x — Xr and A = \J uj 2 — y? . For r > R, 



(2.15) 



<Ptm(r) = ¥>i(r) = B' t hp(y) 



B,r l 



d(r 2 Y 



(2.16) 



where y — ijr and 7 = WfJ^ — uj 2 . Here remember that the condition fio > u> > y is necessary for the existence of 
the stable Q-ball § @. 

First, we consider the 1=1 case. In this case the solution becomes 



(p<(r) = <pi m (r) = -i^-sin(Ar) + Arcos(Ar)) 



and 



tpim(r) = Pi(r) 



El 

2r 2 



7r 



(2.17) 



(2.18) 



The relation between the prefactors, A\ and Bi, is determined by the matching conditions of the field on the wall. 
As the potential Vg is continuous, one should require the ipf (R) = cpf(R) and ^{"{R) = <p'^(R). Each requirements 
imply 



and 



Thus, one finds 



(1 + jR)e^' R B 1 = 2 sin(A-R) - XRcos(XR) A u 

Ax (~2XRcos{XR) + 2sin(Ai?) - (Ai?) 2 sin(Ai?)) = B x (^{iRf + 7-R + l) e 
-2acosa+(2-a 2 )sina (3 2 + 2/3 + 2 



sma — acoso 



1 + /3 



(2.19) 



(2.20) 



(2.21) 



where a = XR and (3 = jR. For Ry = 1 and i?/Lto = 4, the eigenvalue wi is cji ~ 3.6i? _1 . For = 1 and i?/io = 5, 
the eigenvalue wi is u>x — S&R^ 1 . The perturbation is permitted only in the case of the discretized frequency. Note 
that 1 = 1 modes of the surface wave do not exist in the Coleman's study. As I = 1 modes can be easily created than 
1 = 2 modes in general, the existence of the I = 1 modes is important in the context of cosmology. 

Next, we give some conserved quantities which specify the feature of the perturbation. The charge of the pertur- 
bation is 



5Qr 



d xjo = u>x / d xip 1 



2 ui 



R 



+ -a + -cos(2a) + a-sm(2a) + 
2 2 4 V 



sin(a) — acos(a) 

iTp 



-J 



The angular momentum J lm is 



Jx m = J d 3 xT aip = muox J d 3 x(p 2 = mSQ. 



Finally, we consider the 1 = 2 case. From the eqs. (2.15) and (2.16) the solutions are 



ip<(r) = ip 2m (r) 



An r 



3sin(Ar) — 3Arcos(Ar) — (Ar) 2 sin(Ar) 



(2.22) 



(2.23) 



(2.24) 



and 
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p>(r) = ^ 2m (r) = ^(3 + 3 7 r + 7V) e" 

The matching conditions, tpf(R) = ^{R) and ip'f (R) — </?'^(i?) at r = i?, give the relation 

(3 — a 2 )sina — 3acosa 3 + 3/3 + /3 2 

(-9 + 4a 2 )sina + a(9 - a 2 )cosa ~ 9 + 9(3 + 4f3 2 + /3 3 



(2.25) 



(2.26) 



For /j,R — 1 and /j,qR = 4, the solution of ui does not exist. This comes from the fact that the bottom of the potential 
Ve raises as £ increase and then decrease the number of the bound state. For fiR = 1 and /j,qR — 5, one has the 
solution UJ2 — 4.66i2 . 

III. EXCITED STATE OF LOCAL Q-BALLS 

So far, we have analyzed the perturbation on the global Q-balls. In this section, we perform the perturbation 
analysis on the local Q-balls coupled to a U(l) gauge field. This analysis will be significant to cosmology because such 
object might be able to supply a seed of astrophysical magnetic field. 

A. Local Q-balls in the Ground State and Phase Transition 

The basic features of the local Q-ball was studied by Lee et al in Ref. |L4j . We consider the theory with a complex 
scalar field <j> coupled to a U(l) gauge field A^. The Lagrangian density is 



1 



p pw 



(3.1) 



where F^ v — d^A v — d u A^. The conserved current and the energy momentum tensor which induce the conserved 
quantities are 



and 



T, 



1 r 
2 



3n = 2 + ieA^cp* - 4>*{d^ - ieAJ 



(d,j, - ieA^idy + ieA v )(f>* + [d v - ieA^^d^ + ieA^ 



respectively. 

To find the Q-balls in the ground state, we seek for the solution such that 

<P = tp(r)e lut and A^ = A {r)S f , 
Using the above configurations, the Lagrangian becomes 

L = 4tt J drr 2 

where g = uo — eAo(r). The conserved charge associated with the U(l) symmetry becomes 



1 J2 , 1 J2 1 1 „2 .2 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



Assumed the thin wall approximation, as well as the case of global Q-balls, the configuration of the scalar field, ip, is 
roughly expressed by the radius, R, and the field value at the minimum of the potential, er+. f| 



§Rigourously speaking, the value of the scalar field ip for r < R should be determined by the variational principle with 
respect to the expression for the total energy E [Q. One can approximate it as <r+, however, if only the leading-order term is 
considered. The radius R is decided so that the energy of the non-trivial configuration takes minimum value. 
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It is given by 

<p(r) = a+0(R-r). 

In this case, we can solve the equation of motion for g [^4[ . The solution becomes 



e 2 Q 



9(r) 

g(r) = to 



4ttr)^ sinh(e(T + r) 



sinh(e(7 + i?) 
e 2 Q 
Aixr 



for r < R, 



for 



r > R. 

Inserted the above solution into the eq. (3.6), the relation between Q and u> is determined by 

e 2 Q 1 

LO = 

AttR 1 - 

where x = ea + R. The total energy of the Q-ball becomes 



tanha; 



(3.7) 

(3.8) 
(3.9) 

(3.10) 



E = \loQ + 47r / drr 



rl 



U(<p) 



e 2 Q 2 



1 



tanhx\ _1 



■4\/2 



7T(7j 



U-R 2 



Air , 



-U+R A 



8nR V x J 

If ea+R > 1, as will be discussed later, the energy can be approximately written by 

e 2 Q 2 



E 



8nR 



-4\/27i 



rj All o 

U-R 2 + —U+R 3 . 



(3.11) 



(3.12) 



Let us consider the phase transition. We note that our study is different from the previous one by Ellis et al p5| . 
In Ref. jn|, they discussed the phase transition which occurs above the critical temperature, and took no account of 
the existence of the critical charge. We discuss the phase transition which takes place below the critical temperature 
and above the temperature at which the standard nucleation [[Toll starts. 

As the universe expands, the temperature of the universe cools down. At the critical temperature T c , the potential 
has the two degenerate minima at ip = a + and ip = ; U+ = U(Q) = 0. In this case, the total energy becomes 



E 



e 2 Q 2 
8irR 



aV2tt<t + ^/U~R 



(3.13) 



There exists the upper limit on charge and radius of Q-ball 
by 



R* = 



e 2 Q 2 



1/3 



v64v^7r 2 o- +v ^ 

Inserted this equation into the eq. (3.13), the minimum energy becomes 

E* = E(R, 



From dE /dR\u = R t = 0, the radius of Q-ball is given 

(3.14) 

(3.15) 



2 . 25/6^1/3 

Here note that E* is the monotonically increasing function of Q. To construct stable Q-balls, however, the total 
energy E* cannot exceed the energy for Q's free particles with the mass /i, otherwise the Q-balls will decay into Q's 
free particles. Thus, from the condition fi = 9£ l */9Q|g = Q max , one can obtain the upper bound on the total charge 



"Although the configuration of the scalar field is approximately expressed by the step function, we considered the surface 
term with the finite width here. Readers might feel a confusion. However, after all, if we assume ea+R > 1, one see that the 
inconsistency of the approximation can vanish. 



G 



0max = 7^^=- (3 ' 16) 



The upper bound on the radius -R max of Q-ball is also given by 

o _ 1 f 



4V2e 2 a+JTI. 



(3.17) 

We now consider the temperature T q (< T c ). We assumed that T q is higher than the temperature at which the 
bubble nuclcation starts. In Ref. it was pointed out that the Q-ball has the critical charge Q c . If the charge 
exceeds this critical value, Q-ball can expand. If the charge is smaller than the critical one, there exists the range of 
the total energy in which the Q-ball is bounded. At T = T q , if Q-ball is able to have the critical charge, it is possible 
that the phase transition is precipitated by the Q-ball growing up into the macroscopic size. For the Q-ball with 
critical charge, various quantities are decided by the variational principle so as to satisfy the conditions 

dE/dR\ R = RatQ=Qc = d 2 E/dR 2 \ R=R ^ Q=Qc = 0. (3.18) 

Resultant expressions for Q c , R c , and E c are respectively obtained as 



(3.19) 



V2 \U+ 

and 



3 a+JU- 

Rc = ^= + J T | , (3.20) 



E c = 18V2n + T. ■ (3.21) 
\U+r 



As in the case of T = T c , the energy must satisfy 



r- / U- N 1 / 2 

E c /Q c = V3ea+(—) < fi. (3.22) 



According to Kusenko |13|, the large charge can be attained soon and then the almost of Q-balls has the maximum 
charge Q max - If Qmax > Qc, Q-ball can expand as soon as the temperature becomes lower than the critical one 
because the growth rate of the charge is ~ r] v n J (4 : TrR 2 )y / T/m ip ~ 10 _7 GeV S> H ~ 10 _15 GeV for TeV scale. Thus, 
the expansion of the Q-ball will start before the nucleation of the critical bubble. To compare Q c with <5 max , let us 
consider the ratio 



= 12x/3e 3 r±f^) 3 , (3.23) 



^\u + \ 1/2 

where we used the approximation, U^(T q ) ~ f7_(T c ). The condition Q c /Q ma x < 1 implies 

2^V3ea + (^- l ) 1/2 < M . (3.24) 

Comparing eq. (3.22) with eq. (3.24), one can see that the most strict constraint is the latter one (eq. (3.24)). 

Here we must remember that ea+R > 1 was assumed in the derivation of the eq. (3.12). The conditions e<7+i? max > 
1 and eo+Rc > 1, therefore, must be satisfied, so that the coupling constant e is given by 

To see the details, we consider a potential 

1%) = i M V - l -M^ + (3.26) 
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where fi, M, and A depend on the temperature of the universe. At the near critical temperature, one can take 
y = A 1 / 2 (/i/M) = y c + e, where e is a small non-dimensional quantity, and y c — l/\/3 is the value of y at the critical 

temperature, decided by U+ = 0. The leading order of these quantities become U+ ~ a "+ A ^ y " U- ~ M~T~' ° + — 
^|^,cr_ ~ -y-, and jU ~ Vc^- Substituting the above equations into eq. (3.24), one can obtain the constraint on the 
coupling constant e; 



(3.27) 



Furthermore, from the eq. (3.24), the range where the present approximation is valid becomes 



1 A^| e | < e < -^A 1 ^. (3.28) 

We can see that the constraint (eq. (3.27)) can work for a sufficiently small value of e. 

We investigated the phase transition caused by not global Q-ball but local Q-ball. To accomplish the phase transition 
caused by local Q-ball, we showed that the coupling constant e must be too small. In other words, if the model satisfies 

3 3 / 2 e 2 

- ¥Jr <X\e\, (3.29) 

the charge of Q-ball can exceed the critical value as soon as the universe cools down under the critical temperature. 
In this situation, the phase transition will go as follows. By solitosynthesis, Q-balls with maximum charge Q m ax is 
created above the critical temperature. If the coupling constant satisfies the above constraint (eq.(3.29)), the charge 
of Q-ball exceeds the critical value, Q c . Then all the Q-balls must expand to macroscopic size at the temperature 
lower than the critical one. Such phenomena always happen before bubble nucleation starts. Thus, the new picture 
of the phase transition due to gauged Q-balls could become real one. 

We note that the present scenario of the phase transition is slightly different from the Kusenko's one. In Ref. [fl3|| , 
the accretion of charge which happens only below the critical temperature was considered. However, it is conceivable 
that such mechanism will take place above the critical temperature too, and we have taken account of it. 

B. Stationary Perturbation on Local Q-balls 

Let us consider the stationary perturbation on the local Q-balls. The perturbed quantities are written by 

4> = (<p + V^Wi^e*"* (3.30) 

and 

A, = (°>i, + (3.31) 

We assume the following ansatz on perturbation of the vector potential: ^'Aq = ^ A r — 0, and ^' Ai = 
: ^ge m (r)d I (Y em - V^nO where suffix I runs over 9, (p. 
The Lagrangian for the perturbations becomes 



<)L = I d 3 

d :i x 

Thus, the equations of the perturbation become 



I / / \2 -^22 1 TTtt ( \ 2 ^(^ 2 



/ I \2 2 2 2 i 1 *(* + 1) /o oo^i 

(9lm) ~ C <Po9tm + e 9im<P0<Pim n ' ( 3 ' 32 ) 



9im + e 2 (plge m = ^eLp fe m (3.33) 



^Note that all other types of the expression cannot couple with 
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and 



d 2 2 d 2 TT „, . 



+ 1) 



ei(t+l) 
<pe m — n fogem, 



respectively. 

The total charge, the angular momentum and the total energy are given by 



Q = J d 3 xj = J d 3 xg (<fl + <p 2 hn ) = Q + SQ r 



(3.34) 



(3.35) 



d xTq v = ml d xgo 



(3.36) 



and 



E 



= \ojQ Q + j d 3 



1 i 1 If 

2<Po + U(<po)\ + -ujSQ m + - d 3 xg 2 (p 2 n 



=: E + SE, 



(3.37) 



where SE = (l/2)ujSQ m + (1/2) / d 3 xg 2 Lp\ m . 

As we stated in Introduction, the perturbed Q-ball has the magnetic field. The averaged magnetic field B l = 
e l i k djAk over the angular direction is 



B 2 



dnB 2 = 2£(£+l) 



g' 



(3.38) 



We consider only the t = 1 mode of excited Q-balls, since it is unrealistic that the Q-balls with much higher excitations 
actually affect. For the order estimation, we assume ecr i? <C 1. In this case g$ <~ u — e 2 Qo/ATtR. Furthermore, we 
also suppose the small coupling constant, e <C 1, and expand the part of perturbation in powers of e, 



(°), 



¥?i + e (1 Vi and g x = g tn 



(0) .gi+e (1) 5i. 



The equations in the order O(e ) become 
d 2 2d 



The solutions of these equations become 



m u>? = Ax- — - and ( "Vf = B x r- 

dr r 



dr 2 r 



and 



( %i = 



Assumed that a< 1 and (3 -C 1, the expression of 8Qx becomes 

5Qx = ^A 2 ^. 
Up to the order of 0(e), one can obtain the relation 

~ ^ i? 2( °Vi- 

The energy induced by the magnetic field becomes 

[R 

E M ~ 4tt / drr 2 B 2 (r) 
Jo 
16tt 



(3.39) 
(3.40) 

(3.41) 
(3.42) 
(3.43) 
(3.44) 



-e 2 a + R 3 {^ Vl (R)) 2 



16tt 



e 2 A 2 a 2 + R 5 X 6 



407T 2 a 2 + 



(3.45) 
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where we used 5Q ~ Qq. The mean magnetic held becomes 



B : J2E m /^-R3 ~ eQi/^A) 172 , (3.46) 



3 ~ ^ u 

where p = luR. For the potential (3.26), the order of the magnitude becomes 



A 



IV. AN IMPLICATION INTO COSMOLOGY - PRIMORDIAL MAGNETIC FIELD - 

In this paper, we investigated the perturbation on the ground state of Q-balls. Under the thin-wall approximation, 
we solved the equation for the part of perturbation analytically in the global Q-balls. For the local Q-balls, we 
solved them by expanding in powers of coupling constant e. Furthermore, we estimated the mean magnetic held for 
each Q-balls. In the context of cosmology, we need to average it. We assume that gauged Q-balls can be copiously 
produced. Since it is natural to suppose that the excited states of the Q-ball obeys the thermal distribution in such 
situation || , the cosmological mean magnetic field generated by the excitation of the local Q-ball can be evaluated 
by multiplying the factor e~^ 1 ^ 2 ^ SE , where SE is the energy of the excitation and is given by 

1 Ana 
5E ~ 2 x -uSQx ~ loQo ffi. (4.1) 



2 



Thus, the mean value becomes 



where £ := (l/2)f35E = 



At the endpoint of the phase transition, the magnetic field becomes 



fimax — U "J 72~3~7T72 ™2 e 



(B)n f ~ (^p) (B) 



R f ) n " / 6 2 e 3 p 1/2 %i 

\ \( fb \- 2 / e \-3/ p n-i/2/ T- 



~l"xf-)f-)fA)1-)'°W( lf )W-800*mf— Gauss. (4.3) 
VlOO/ VI. 0/ V10- 3 / Vo.l/ U.0/ VIOOGeV/ P L U.0A0.1/ tJ v ' 

where Rf denotes the size of Q-ball at that time and /& = RfHf. Since the Reynolds number becomes Re ~ 10 12 
at the endpoint of the phase transition, turbulence for the electroweak plasma occurs over the scale of the radius 
of domain and then the energy of the magnetic field is equipartitioned with the energy of the fluidp|. Thus the final 
field strength becomes 

B(R f ) ~ p^v ~ gl /2 vT] ~ l0 24 ( I ^- v ) 2 Gauss, (4.4) 

where p and v are the density and the velocity of the fluid, respectively. 

If the coherent scale is comoving, the number of magnetic domains inside the galaxy scale is N ~ 
10 10 (10~ 3 //b)(T//100GeV). Thus, averaging over the galaxy scale JTsfl, one can obtain the present mean value 



-2 



(-E> now ) galaxy ~ AT3/2\„ ) ^C^/) 



4 now ' 



"Although the amplification due to the turbulence is not clear at present, it has been used in the generation mechanism of 
the primordial magnetic field during the first order phase transition j?j [^| . 
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The above value satisfies the onset condition of the dynamo mechanism. 

In the above derivation, we assumed for simplicity that the evolution of the magnetic field goes only by the red-shift 
effect. Although the definite scenario of the evolution has not been composed yet, some studies related with this 
problem have been done. Finally, we will refer to the studies which deal with the details of the evolution of the 
magnetic field. In Rcf. JiT| , the damping effect due to photon diffusion around the recombination era were studied. 
In Ref. (^0|, Olesen et al showed that an inverse cascade happens if the power of the initial spectrum is larger than 
—3, and this means that the coherent scale of the magnetic field will be able to extend. 
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